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The solid-on-solid model provides a commonly used framework for the 
description of surfaces. In the last years it has been extended in order 
to investigate the effect of defects in the bulk on the roughness of the 
surface. The determination of the ground state of this model leads to 
■ a combinatorial problem, which is reduced to an uncapacitated, convex 

C minimum-circulation problem. We will show that the successive shortest 

path algorithm solves the problem in polynomial time. 

ON 



1 Introduction 



S3 . 

In the past twenty years there has been a fruitful collaboration between theo- 
retical physicists and researchers from the field of combinatorial optimization 
on the determination of the global energy-minimum (ground state) in various 
^ ■ spin-glass models. The observation, that the ground state of the Ising-model 

can be found by solving an certain maxcut-problem, has drawn the attention 
of physicists to the methods developed in graph theory . Conversely this ap- 
plication has animated the research on the maxcut-problem |3]|| itself. It the 
same spirit it was shown, that in the context of the Random-Field-Ising-model 
the maxflow-problem arises ||. 

The solid-on-solid (SOS) model, on which we will comment here, is a micro- 
scopic approach to describe the aggregation of atoms on a crystalline sub- 
strate. In contrast to the above mentioned Ising-models, where the dynamic 
variables can take on only two values, in the SOS-model these variables are 
integer valued. At first sight this fact seems to complicate the optimization 
problem. However, similar to the Ising-model, the underlying structure can be 
represented by a graph. We will show in this paper, that in the planar case 
the problem of minimizing the surface-energy can be reduced to the problem 
of finding a network-circulation, that minimizes a convex objective function. 
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The growth of surfaces is a subject of intense research in chemistry and solid 
state physics. A deeper understanding of the relevant processes has applica- 
tions in the fabrication of semiconductors and various other fields where the 
precise shape of the surface is crucial ||. In the simple approach taken by 
the SOS-model the lattice structure of real crystals is taken into account by 
subdividing the surface into discrete sites (e.g. unit-squares) on which par- 
ticles, modeled by unit cubes, can aggregate [BUT]. There are several possible 
parametrizations of the model. The one which is commonly used for theoretical 
investigations yields the surface-energy 

H(h) = - h 3 )\ (1) 



Here hi G Z is the height of the surface at site i and (ij) indicates all pairs 
of nearest neighbor sites on e.g. a square lattice, mimicking the crystalline 
structure. Obviously the ground state of this model is flat, i.e. hi = const for all 
i. Thermal fluctuations lead to a roughening transition at some temperature 
above which the height-height correlation function increases logarithmically 
with the distance. 



One assumption of this model is the flatness of the underlying substrate, 
i.e. a bulk-ordered crystal. Due to the presence of quenched defects one should 
consider the intriguing question of what the effect of a disordered substrate 
on the thermodynamic (and non-equilibrium growth) properties of the surface 
might be |14}1 . It has been suggested |y| to introduce a disordered substrate 
into the model (|I]) by providing each of the height variables with a real valued 
offset di G R, which is generally chosen to be from [—0.5, 0.5]: 

hi := rii + di with rii G Z. (2) 



This model has some new features not contained in the original one. The in- 
fluence of the disorder tends to roughen the surface also at low temperatures. 
In particular a new phase has been predicted |L5,iI|, where height-height cor- 



relations diverge more strongly with the distance than logarithmically. Since 
this is a low temperature phase the properties of the ground state of (H) with 
(0), which is no longer trivial due to the competition between the elastic term 
(hi — hj) 2 and the random offsets di, are of crucial interest. 

The paper is organized as follows: Section introduces the reduction to a 
convex network-flow problem and in section |3] we will then show that (p]) can 
be minimized in strictly polynomial time using the techniques introduced in 



1 1|. We use standard notation from graph theory, and refer the reader to 



for the definitions. 
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2 Reduction to a network-flow problem 

Let D = (V, E) be a directed graph. The optimization problem we consider is 

min H(n) = J2 d d i + n i) ~ ( d i + n o)f ( 3 ) 
(ij)eE 

s.t. 

rii G Z forall i G V. 

Using the definitions dij := di — dj and := — the objective-function 
can be written as 

H(X)= (dn + Xijf. (4) 

We intend to reformulate the original optimization problem in terms of the 
variables x^j. Obviously we must demand x^ G Z for all edges (ij) G E. 
Furthermore, since the Xij describe potential-differences in the scalar field 
given by a set of height variables rij , it is clear that the sum of the x^ along any 
oriented cycle on the surface has to be zero. The following theorem establishes 
this constraint formally in the special case, where D is a planar graph. 

Theorem 1 Let D = (V, E) be a planar, directed graph and D* = (V*, E) be 
the associated dual graph. There exists a set of node-variables ni, ... raiyi such 
that X = (xij) E = (ni — nj) E if and only if X satisfies 

x^ — x ji = f or Q tt dual nodes i G V*. (5) 
(ij)eE (ji)eE 

A proof of this theorem can be found in Q Chapter 5. 
Given a vector X and a node i we call the number 

9i( x ) '■= x ij ~ X H 
{ij)&E (ji)eE 

the excess of a node i. A complete formulation, which, for planar graphs, is 
equivalent to (|3|) is thus given by 

wmH(X) = £ (dy + xy) 2 , (6) 
(ij)eE 

s.t. 
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Xij G Z for all edges (ij) G E, 
gi(X) = for all dual nodes % G V* . 



(7) 
(8) 



A vector X satisfying these equations is called a circulation The problem 
is a convex network-flow problem in the dual graph without any capacity 
constraints (the flow-variables on each edge can take on arbitrary positive 
or negative values). We note that this reduction, as well as the results in the 
following section, more generally hold for any objective function, which can 
be written as 

h(x) = J2 

(ij)EE 



where the On are convex functions with a minimum. 



3 Algorithm and complexity 



The algorithms which solve the convex flow problem [J IJ are not strictly poly- 
nomial if there are no capacity constraints. We will show in this section that in 
the special case of problem (|6[) the time complexity of the successive shortest 
path algorithm for convex flow can be bounded by a strict polynomial. For 
this purpose we start with a brief review of the algorithm, a more detailed 
description can be found in 0. 

The underlying philosophy is to search the solution in a subset of Z E , in 
which each vector X yields a lower bound for the optimum of (]£]) but does 
not necessarily satisfy the flow-constraints (|8|). The flow balance at all nodes 
has to be approached iteratively by the algorithm. To characterize this subset, 
one introduces the notion of dual feasibility: 

Definition 2 Let H(X) = J2(ij)eE @iji x ij) be a sum of functions of the 

components of X. A vector X is called dual feasible, if there is a number 
Hi for every node (potential-function) , such that for all edges (ij) G E the 
reduced cost inequalities 



C lj{ X ij) '■= 9ij( X ij + 1) - ®ij( X ij) - + TTi > (9) 
Cjii^j) '■= 9ij( X ij - !) _ ®ij( X ij) + - TTi > 

are satisfied. 

Theorem 3 Let X be a dual feasible vector and let X opt be an optimal solution 
of problem (g), then H(X) < H(X opt ). 
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The proof is given in |T| chapter 9 and chapter 14.3. An immediate consequence 
of this theorem is 

Corollary 4 A dual feasible vector X satisfying Q), (|<|) is an optimal solution 
of problem (Qj. 

In order to apply the successive shortest path algorithm we determine a dual 
feasible initial solution. Here we can choose the vector X s 6 Z* E , whose compo- 
nents minimize the functions 9ij(x) and the potential function n := (7Tj)l^ = 0. 
Note that, when tt = 0, the reduced costs (|9]) simply measure the change of 
the objective function triggered by a change of the (zj)-component by plus 
resp. minus one. These are clearly all positive if X minizes every summand of 

©• 

The initial solution does in general not satisfy the flow-equations (H). The 
excess g%{X) of all nodes is then balanced by iteratively augmenting X along 
paths between supply-nodes (gi(X) > 0) and demand-nodes (gi(X) < 0), 
which are shortest paths with respect to the reduced costs. The variables Xij 
associated with the edges along such a path are changed by exactly one unit 
(the sign depends on the orientation of the edge in the path). 

Due to the non-linearity of the objective every update of X leads to a new set 
of arc-weights. However, the convexity of the individual functions 9ij ensures, 
that the inequalities (§) still hold after every augmentation, when c^^x^) is 
replaced by cjAxij + a) with a G { — 1,0,1}. Therefore the algorithm works 
correctly because it always maintains dual feasibility of X with respect to the 
current set of arc-weights. 

procedure successive shortest path 
begin 

X= min{H(X)\X eZ E } 

7T = 

while there is a node s with g s {X) > 
compute the reduced costs c*(X) 

determine the shortest path distances d(i) from s to all 

other nodes with respect to the reduced costs 
choose a node t with gt{X) < 

let w(s,t) denote the shortest path vector from s to t in D 
X = X + w(s,t) 

IT = TV — d 

end 
end 

The complexity of the algorithm is 0(G(X s ))S(n, m). Here S(n, m) is the time 
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needed for a shortest path search in a network with positive arc-weights and 

G(X S ) :=^EI#(*.)I 



counts the number of paths needed to balance the excess of all nodes. The 
following lemma then yields an upper bound for the value of G(X S ). 

Lemma 5 The vector X s e 1> E that minimizes H(X) (eq. ^) satisfies G(X S ) < 
\E\. 



PROOF. If the integrality- as well as the flow-constraints are neglected the 
global minimum of H(X) = J2(ij)eE(dij + Xij) 2 is attained at —d := (—dij) E 
with H(—d) = 0. If we now take into account the integrality-constraints it 
follows from the convexity of the (dij + Xij) 2 , that the minimum xfj G Z is 
attained in either ["— dy] or |_— dyj ■ We thus find 

I4-K-)I< L 

From theorem [I] we obtain G(—d) = and thus 

G(X S ) = G(X S ) - G(-d) < X; K-{-dij)\ < \E\. 

(ij)eE 



4 Computer experiments 



We implemented the successive shortest path algorithm for the convex-flow 
problem. The program needs less than 80 seconds on a Sparc Workstation 
to solve instances of grid-graphs with 128 x 128 nodes and randomly chosen 
offsets to optimality. Grid-graphs of size 200 x 200 nodes take about 10-15 min- 
utes. We used C++ for our implementation and made use of the LEDA-class 
library [12]. The physical results, which have been obtained with this program, 
together with their theoretical interpretation will be reported elsewhere [O. 
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